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ABSTRACT 


The acceleration o£ energetic ions in interplanetary magnefcosonic fast** 
inode shock waves is studied via analytical modeling and numerical simulations* 
An analytical model that combines both the shock drift and compress Iona l 
acceleration mechanisms is presented. The analytical predictions of the model 
are shown to be in good agreement with numerical simulation results. 



I. INTRODUCTION 


Observations during tins past two decades at 1 All and in deep space have 

established a causal relationship between the passage of interplanetary shock 

waves and the large enhancements of energetic proton intensity that are frequently 

observed around the time of shock passage [Snrris and Van Alien, 1974; Peases 

et al . , 1979, and references in both papers). 

Axford and Reed [1963) first suggested that the above relationship is due 

to protons being directly accelerated in interplanetary shock waves. Since 

then several types of interplanetary shock acceleration mechanisms have been 

discussed: Compression between the shock front and upstream magnetic field 

irregularities [Fisk, 1971]; A gradient | B { drift at the shock front in live 
> *► 

V x B electric field in the shock rest frame [Chen and Armstrong, 1972, Sards 
and Van Allen, 1974; Armstrong et al . , 1977] ; and compression between upstream 
and downstream magnetic field irregularities [Fisk and Leo, 1980). 

The purpose of this paper, the first of three companion papers, is: 

1) to construct a model which combines the above acceleration processes and 

2) to present analytical expressions for energy changes, final pitch angles, 
and acceleration times for particles transmitted and reflected by oblique and 
perpendicular magnetosonic fast-mode shocks. The model presented here is used 
in Paper XI to calculate the differential energy spectra of ions accelerated in 
corotating Interaction region (CT.R) shock waves; and compared In Paper I'll to 
high time resolution observations of CIR accelerated protons. 



II. SHOCK GEOMETRY 


Observations of energetic (~ i MeV) proton events associated with shock 


waves in deep space {Pesses et al ., 1979J and at i All {Snrris, 1973, and refe- 


rences therein] show that these shocks are apparently all of the magnetosonic 


fast-mode variety. No observations of energetic proton events associated with 


either magnetosonic slow-mode or Alfven shock waves have, to my knowledge, been 


reported. Therefore, the acceleration model to be developed here will be Cor 


only the fast-raode shocks. 

>■ 

The magnetosonic fast-mode shocks to be considered move with a velocity V s 
with respect to the upstream plasma rest frame. The plasma mass density up- 


stream of the shock is p\ and downstream pg. The shocks are planar in the y - z 
plane of Figure 1, The shock front unit normal vector n is directed along the 

positive x-axis, and ( 1 ^ 2 ) ** 8 fc * ie acute angle between *n and Bi (B;>) * the up- 

stream (downstream) magnetic field vector. The. hydrodynamic shock strength is 
P 2 P 1~ ' and t,ie magnetic shock strength is | B 2 | | I” 1, (<* N) • The 
motion of the magnetized solar wind in the shock rest frame results in a Vj. x 8 

electric field. In Figure 1 this E is parallel to negative y-axis. From the 

■>* 

continuity of the tangential component of E across the shock front the upstream 


(E^y) and downstream (Egy) electric field-vectors are equal. In the model it 
is assumed that V 8 , n, B^ and do not vary with time or space. 


III. SUGCK-ASSOCIATED ACCELERATION PROCESSES 
A. Classification 

A general classification scheme for charged particle acceleration mechanisms 
has been developed by Northrop [1961)* He shows that it a given reference 
frame the time (t) rate of change of a non-relatlvistic charged particle’s 
kinetic energy (T) averaged over one particle gyroperiod, <dT/dt>m, is 
given by Equation (1) 

I It 

<££>« « qG*<E>... + M it JJlL d* 0(m 2 /q^> (1) 

dt w fit 

-K 

In Equation (1) q is the particle charge, m is the particle mass, G is the 

+ 

particle guiding center velocity, <E>m is the average value, over one gyroperiod, 

4 ¥ 

of the electric intensity vector E at the position of the guiding center, B is 

the magnetic inductance vector at the position of the guiding center, p is the 

magnetic moment of the particle (evaluated In the guiding center rest frame), 

and 0 means on the order of magnitude of. Term I in Equation (1) is the time 

rate of change of the particle energy due to work done by the electric field on 

> 

the guiding cent****. Term II is an induction effect of a time dependent B and is 

> 

the time rate of change of the particle energy due to the curl of E acting 
about the circular particle gyro orbit. 

Shock drift acceleration is included in Term I, and compressions! acceler- 
ation in Term II. Northrop [1963] notes that both betatron acceleration and 
the type of acceleration first discussed by Fermi [19A9] are included In Term IT. 

B. Reference Frame Dependence 

Interplanetary shock acceleration models calculate the particle energy 

gains in various frames: the shock rest frame, the upstream or downstream 
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plasma rest trnmaa, and tha bull electric field frame where both the V 8 x B and 
♦* 

an/Bt elactric fields are aero* 

1. Shock Rest Frame 

4 > 

In the frame la which the shock ts at rest the V)b| drift motion and also 

» ► 

drifts due to changes across the shock front of the direction of B and the R x 
> 

B drift velocity drift velocity result in a net displacement of the particle 

* > 

guiding canter in the V s x B electric field potential. This type of shock 

>• > 

acceleration is called shock drift acceleration (previously called V x 11 
acceleration). 

► 

As there is no SB/3t due to the shock’s motion in this frame, there is no 
"induction" acceleration from the shock front. However, charged particles are 

*h 

also accelerated in this frame from the DB/Bt. curl of 12 produced by moving 
magnetic field irregularities. Particles baekaentfeered toward the shock by 
approaching upstream irregularities gain energy while particles bnckscottered 
toward the shock by receding downstream irregularities lose energy. The mag- 
netic field irregularities are converted by the plasma bulk motion so the up- 
stream irregularities are approaching faster than the downstream irregularities 
are receding. This divergence in the velocity of the irregularities at the 
shock front results In a net energy gain due to compression. (Ax ford et aJL. , 

1977 j Bell, 1978) . 

2. Upstream Plasma Rest Frame 

In the upstream plasma rest frame charged particles arc accelerated by the 

curl of R induced by temporal variations in B duo to both the motion of the 

shock front and approaching downstream magnetic field irregularities. A graph ie 

way of looking at the acceleration of the shock trout hi this frame N given 

in Figure 2a. The energy gain Is due to particles gyrovelocity having a component 

~ 6 - 



+ *► 

which Is parallel to the downstream V x B electric field. 

3. Downstream Plasma Rest Frame 

In the downstream plagrap rest frame charged particles are accelerated by 

+ + 
the curl of E Introduced by temporal variations in B due to both the motion of 

the shock front and approaching upstream magnetic field Irregularities. A 

graphic way of looking at the acceleration at the shock in this frame is shown 

in Figure 2b. The energy gain Is due to the particles’ gyroveloclty having a 

4 - > 

component parallel to the V x B electric field in the upstream region. 

4. Null Electric Field Frame 

In the null electric field frame (which moves parallel to the shock front 

^ •> + 

in the shock rest frame with a speed (V s .1itant|ii -|v 8 x n|), both <E>w and 
9B/ 8t are zero at the shock front. Hence in the null electric field frame par- 
ticles are accelerated solely by the electric field produced by moving magnetic 
field irregularities. 

To summarize the above discussion, there are two basic energetic particle 
acceleration processes associated with fast-mode shocks: the shock drift mecha- 
nism which is present in the shock rest frame, and the "inductance" mechanism 
which at Che shock front is present in both the upstream and downstream plasma 
rest frames, and which at moving magnetic field irregularities (compression) 
mechanism which is present in all four frames. 

Several models of interplanetary shock acceleration have neglected one 
of the above acceleration processes. Fisk's [1971) snow plow model assumes 
that particles are accelerated by compression between the shock front and up- 
stream magnetic field irregularities. This model explicitly assumes that 
reflected particles do not undergo shock drift acceleration, and implicitly 
assumes that transmitted particles gain no energy. 
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Both the energetic storm particle event models of Scholar and Morfiil 
[1975] and the corotating particle event modal of Palmer and Gosling [1978] 
assume that particles are accelerated by multiple reflections off the shock 
front. Both these models assume that reflected particles are shock drift accele- 
rated but that transmitted particles undergo no energy gain. 

The shock acceleration models of Sarris and Van Allen [1974] , Armstrong 
et al . [1977] and Pecker [1981] assume that in the shock rest frame both ref- 
lected and transmitted particles are shock drift accelerated. These models 
are basically concerned with quasi perpendicular (tl ~ 90°) shocks and 
ignore compresslonal acceleration which is not significant compared to shock 
drift acceleration at t|i ~ 90° » 

The corotating particle event model of Fisk and Lee [1980] does not negate 
either shock drift or compression. Their model assumes that particles. In the 
null electric field frame* are accelerated by compression, and thus in the 
shock rest frame particles are. accelerated by both the shock drift and compression 


mechanisms. 



IV. SINGLE SHOCK ENCOUNTER ANALYTICAL EXPRESSIONS 


A. Post Encounter Energies 
1. Choice of Coordinate Frames 

The energy a particle gains from a single reflection or transmission from 
a shock can be calculated in any particular frame. However, to calculate the 
total kinetic energy a particle gains from multiple shock encounters it is 
necessary to include energy changes that result from the particle’s being 
backscattered to the shock by moving magnetic field irregularities. If the 
expression for the particle energy gain from a single shock encounter is given 
in terms of the pre-encounter and post-encounter energies in the particle prc- 
and post-encounter plasma rest frame, then the particle backscattering , which 
is assumed elastic in a plasma rest frame, will result in no additional energy 
gains. The backscattering, congressional energy gains come from the relative 
velocity of the upstream and downstream plasma rest frames. 

The energy gain per shock encounter (one encounter is composed of several 
shock crossings) is calculated using the following procedure; First, the par- 
ticle velocity is transformed from the pre-encounter plasma rest frame of the 
particle to the null electric field frame; second, an algorithm is used to de- 
termine If the particle is reflected or transmitted by the shock and what the' 
post-encounter pitch angle is; third, the post-encounter particle velocity is 
transformed into the post-encounter plasma rest frame of the particle. 

It will be shown in the section on numerical simulations that in the null 
electric field frame the pre-encounter and post-encounter value of u when 
averaged over gyrophase are equal. 
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In calculating the energy gains the following assumptions are made; 

•k 

(1) contributions to the B field due to the interacting particles are negligible; 

(2) energy losses due to electromagnetic radiation are Ignorable; 

(3) no collisions occur between the Interacting particles; 

(4) V a , H, N, and ^ do not change during the particles 1 shock Interaction, and 

(5) the ions to be accelerated are well above the thermal distribution. 

Previous analytical studies of the interaction of charged particles with 
interplanetary shock waves have been carried out by several authors for perpen- 
dicular shocks [Shabanskli, 1962; Schatzraan, 1963; Pesses, 1981] and for 
oblique shocks [Hudson, 1965; Alekseyev and Kropotkin, 1970; Singer and 
Montgomery, 1971; Sarris and Van Allen, 1974; Vansl'yed et al . , 1978]. 

In the upstream (downstream) plasma rest frame the null electric field 
frame must move parallel to B^ (B 2 ), otherwise V x E The null electric 

field frame must also move along the field lines at the rate at which the shock 

f 

does, otherwise 38/3t £ 0. Hence, the null electric field frame moves with 
respect to either the upstream or downstream plasma rest frame parallel to the 
magnetic field vector with a speed equal to the projection of the shock velocity 
along the magnetic field lines. 

The transformation from the ipstream plasma rest frame to the null electric 
field frame is carried out by moving with a velocity V^, where 

V 1 * v lx sec ’hA (2) » 

Vj x = V s .n, and a Bi/|Bi|. The transformation from the downstream plasma 
rest frame to the null electric field frame is carried out by moving with a 

•f 

velocity V 2 , where 
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2. Reflection and Transmission Considerations 
The algorithm that determines if a particle is transmitted or reflected in 
the null electric field frame is derived by demanding that the pre- and post- 
encounter particle kinetic energy and angular momentum about the guiding center 
be equal. The particle angular momentum h is given by the cross product of 
the particle gration velocity and gyroradius, 


sign(q) m 2 V, sin 2 a. 


t 


+ 

B 


(4) 


where V^( a, t ) is the particle speed (pitch angle) in the plasma rest frame of 
the particle. When the particle is directed towards the upstream region 0 < 

< 90°. Numerical simulations to be discussed later show that a particle will 
not be reflected if it can conserve angular momentum and energy by being trans- 
mitted. If transmission would violate a conserved quantity a particle would 
be reflected. The condition under which incident upstream particles are 
transmitted downstream by the shock, and incident downstream particles are 
transmitted upstream by the shocks, are given by Equations (5) and (6), respectively. 

1/2 m NV 2 sin 2 a < 1/2 m ( I V, I + V 2 - 2 1 V-, | V cos a ) (5) 

1 i 1 Ll i i i 

1/2 m jrV sin 2 a < 1/2 m ( | V 2 | + V 2 - 2 1 V 2 |v cos a ) (6) 

it i ll 

For particles initially upstream the pitch angle boundary a l>2 between 
reflected and transmitted particles for given values of V*, and N are found 


by equalling the right and left hand sides of (5) and solving for a. This gives: 




( 7 ), 

<«), 


where R * |Vt|V“ 1 . Upstream particles with cos a, > R cannot be overtaken by 

1 * 1 ^ «** 

the shock and so do not Interact with it. Upstream particles are reflected when 
either cos”^R < cti < ctj and R < 1 or when < ot^ < c *2 and R > 1. Upstream 
particles are transmitted when ctj > 03 and also when provided R > 1* 

When R » JT, * a£ and all upstream particles are transmitted. When 
R > /T, the right hand sides of (7) and ( 8 ) are complex numbers; and numerical 
simulations show that all upstream particles are transmitted [Chen, 1975]. 

for particles that are initially downstream Parker [1963] has shown that 
they are either transmitted upstream or do not interact with the shock. No 
initially downstream particles are reflected. Downstream particles with cos 
< JR, where J * KH“1 , cannot overtake the shock and so do not interact with it. 
Downstream particles with cos aj_ > JR aie transmitted upstream. 

The upstream noninteraction n u , reflection r u , and transmission t u coef- 
ficients of the shock as a function of R and N are given in Equations (9a-c) 
These equations are for an initial pitch angle distribution that is isotropic 
with respect to B in the upstream plasma rest frame. 
f 1 cos"”^ R R < 1 

n u * (9a) 

0 R > 1 

tt~*-(cu *•* it n„) R < 1 

r u - . (9b) 

ir"’ 1 (a 2 - a x ) R > 1 

t u - 1 - r u - n u (9c) 


Some upstream particles will always be reflected by the shock provided 
R < /r and N > 1. The downstream transmission and noninteraction i\j 
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coefficients of the shock as a function of JR and N are given in Equations 

(lOa-b). These equations are for an initial pitch angle distribution that is 

*► 

isotropic with respect to B in the downstream plasma rest frame. 

Tj * Tr”^cos”*(JR) JR < 1 (10a) 

- 0 JR > 1 

fld " 1 ” t d (10b) 


3 . Energy Change 

Using the method discussed previously in section IV-A the fractional 
kinetic energy change, (Tfi na i - ^initial) ^-initial » P er shock encounter in the 
plasma rest frame(s) for upstream particles reflected upstream ATp/Tj, upstream 
particles transmitted downstream AId/T^, and downstream particles transmitted 
upstream ATg/Ti, are presented in Equations (11), (12), and (13) respectively. 

£Tb_- 4r(R - X ) (11) 

T i 

AT P - r{r(J 2 + 1) ■ 2 x ‘ 2 J[1 + R(R ' 2 x) - N s 2] 1/2 } (U) 

Ti 

Ml - R{R(J 2 + 1) - 2 J X +2[l + JR(JR-2 x)'S 2 N -l] 1/2 } (13) 

Ti 

where x a c08 <*i and s “ si n Note that in Equations (11-13) AT/T^ dot's not 

+ 00 as >90° (i.e. as R + °°) . All upstream particles are transmitted when sec ^ 

> V V- 1 vlT. Hence, AT R /T> does not approach °°. Expanding the radical in Equation 
i lx * 1 

(12) and taking the limit •> 90° gives 

AT 

l im <*1 + 90 °) ’ (N “ 1)S2 (U) 

and hence ATp/T^ does not approach 00 . No downstream particles ate transmitted 

upstream when sec t[u > V V*"^ J“ . So AT, t /T 4 does not + 00 . 

1 i lx u A 
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Two examples of how shock drift and congressional energy gains are both 
included in the above equation are given below. Consider a particle that is 
reflected back and forth between the shock front end upstream magnetic field 
irregularities. The particle has an initial energy in the upstream plasma 
rest frame. The post-reflection energy in the upstream plasma rest frame Tj 
can be calculated form Equation (11). After backacattering towards the shock 
particle's energy is still T 2 « The energy in the upstream plasma rest frame 
after a second shock front reflection can be calculated from equation ( 11 ) 
again, and so on. Now consider a particle that is transmitted back and forth 
across the shock front. The particle has an initial energy in the upstream 
plasma rest frame* The particle’s energy in the downstream plasma rest 
frame after it has been transmitted downstream (T 2 ) can be calculated from 
Equation (12). After backscattering toward the shock the particle's energy in 
the downstream plasma rest frame is still T 2 • The particle's energy in the 
upstream plasma rest frame after being transmitted upstream (T 3 ) can be 
calculated from Equation (13). After again backscattering towards the shock 
of the particle energy in the upstream plasma rest frame, is unchanged, and 
Equation (12) can be used to calculate the particle's energy in downstream 
plasma rest frame after it has been transmitted downstream a second time, and 
so on. 

4. Maximum Values of AT/T^ 

The maximum value of AT/T^ in Equation (11) occurs for «* op anc * 

Equation (13) for cos «i * JR. In Equation (12) the maximum value of AT/Tj 
occurs for * «2 when R < /FT and for 

cos - RN ” 1 |l - [(N - 1)(J 2 N - l )" 1 (l - NR ” 2 )] 1/2 }, 

when R > /IT. The maximum values of AT/T^ as a function of R, N, and J are given 
in Equations (15) - (17). -14- 


max AT R (R,N) - 4R 2 iffrll {l + [~L- 1 1/2 } < l5 > 

T l N N-l IT 

max A !p, (R,N, J) - R 2 (j 2 + 1 - 2 {l - [(N-1)(1L, - 1)] 1/2 }) (16a) 

T-£ n r 

R < /N 

max Mji. (R,N,J) - R 2 {J 2 + 1 - 2 {l + [ (N-1)(NJ 2 ~1)(1 - I)) 1 / 2 } 

Tjl N r2 

R > /N " 

max & (R,N, J) - R 2 (1 - J 2 ) + 2R 1(1 - I)(l - J 2 R 2 )]1/ 2 (17) 

T* N 

The largest value of AT/Tj[ as a function of N and J is obtained by finding 
the value of R which makes the value of the equations for max AT/T 1 (R,N,J) a 
maximum. The largest value of AT/l’i for R < /N occurs for reflected particles 
and is given by Equation (18). 

max A L= 2(N-1) {l + [N(N-l)" 1 ] 1 / 2 ] R < /N~ (18) 

Ti 

The maximum value of AT/T-j^ for R > /N~ occurs for transmitted downstream 
particles and is found by substituting R 2 = 0.5 N[1 + (1 - 4 3) “l/ 2 J into 
Equation (16b), where 

8 » N _1 (N - 1)(NJ 2 - 1)(1 + J 2 - 2N _1 )“ 2 . 

The maximum value of AT/T-^ as a function of N is, for N = 2, 4.83; for 
N * 3, 8.90; and for N - 4, 12.93. 


(16b) 
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B, Boat Shock Pitch Angles 


1. General Considerations 

The equations for AT /T ^ presented in Section IV-A3 give the particle total 
kinetic energy change per reflection or transmission by the shock. To calcu- 
late thu post encounter pitch angle it ia necessary to know both the poat- 
encounter parallel and perpendicular kinetic energy of the particle. Northrop 
[1963] shows that in a given reference frame the time rate of change of a non- 
relativistic particle’s parallel kinetic energy (Tj) and perpendicular kmetie 
energy (Tj.) averaged over one particle gyro period <dT^/dt>m and <<lTj/dt>M, are 
given by Equations (19) and (20) respectively. 


i i ' 

<gl > w - qV 1| E 11 - p V„ IM + niV^ . f|B + 0(m/q) 


dt 


•f 

v 

IE 


dt 


(19) 


<^IvL > w 
dt 


iLlL + p V , r a M 4* p Vj, . V |n| mVj,, . 4* 0(w»^/q^*) / 2 o) 

uu U u ll U ' ' 


->• > > 

In (18) and (19) V u (E y } is the component of G [<E> W ] parallel to B, is the 

* -► a •> 

E x B drift velocity, B is a unit vector in the direction of B, and a is the 

distance along the B lines of force* 

In the shock rest frame temporal changes In T d are due to the magnetic 

. ► . +• * 
mirror term, pV||0|Bj/3s, and the interaction of the ExB drift with the change 

•> 

in direction of B at the shock front. Time variations in Tj_ in the shock rest 

frame are due. to the magnetic mirror force plus the interaction of the E k B 

drift with the gradients in j B { and Vg at the shock front* 

In the plasma rest frames time variations in Tjj an* due to the magnetic 

* 

mirror force plus the interaction of the E x B drift either with temporal and 

> 

spatial variations in the direction of B. Time variations in Tj are due. to the 
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magnetic mirror force and the interaction of the 15 x B drift with temporal and 

* i *1 

spatial variations in Vg, and spatial variations in |B|, 

In the null electric field frame temporal variations in both Tj and Tj, 
are due solely to the magnetic mirror force. 

2. Calculation Via Transformations 

The post encounter pitch angles ore derived using the same series of coor- 
dinate system transformations utilized in the post shock energy calculations. 

The post-encounter pitch angle in the particles post encounter plasma rest frame 
for reflected particles (ajg), transmitted downstream particles ( p) and trans- 
mitted upstream particles («fp) are given in Equations (21), (22), and (23) 
respectively. 


tan 

a ■ __ 

S 




(21) 


eu 2 

R - X 





tan 

Qt « 

fD J 

S 

R - [1 + R(R ■ 

/r 
- 2 

w 

riPsTjT/T 

(22) 

tan 



S N 

-1/2 

(23) 

tt ** — 

fU R 

- [1 + J R (J 

“r” 

2 

x j "j . 


Expanding the radical in (22) and taking the limit as ^ *► 90°, gives 

tan ct£p » /IT tan ot^ . (24) 

The above equations and Equations (6) - (20) are all averaged over a gyro tropic, 
phase angle distribution. The post-encounter pitch angles show that as expected 
from angular momentum considerations, reflected particles gain energy only in 
their parallel component and transmitted particles gain energy predominantly 
in the perpendicular component. Note that in Equation (22) some transmitted 
downstream particles are directed back towards the shock in the downstream 
plasma rest frame, i.e*» 
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otjj -0 < 90*, However, no transmitted downstream particle has a sufficiently large 
parallel speed to overtake the shock. 

Equations (21) ** (23) are kinematically consistent In that upstream par- 
ticles reflected (transmitted) by the shock always outrun the shock (always 
are left behind by the shock). Likewise downstream particles transmitted by the 
shock always outrun the shock. 

Equations f 11) - (24) give the poet-encounter energies and pitch angles in 
the plasma rest frame. Since all observations are made In the spacecraft rest 
frame, these equations must be transformed into the spacecraft rest frame. 

The expression for post-encounter AT/T^ and cif in that frame can be obtained 
by making the following substitutions 

R + R + <V w i * n)coa ^ jV-jJ " 1 

JRUR + (V w2 . «)l°*Jl Iv.r 1 

where V w i (V w £) is the solar wind velocity upstream (downstream) of shook In 
the spacecraft rest frame. 


C. Particle-Shock Interaction Time 

The particle-shock interaction time can be estimated by calculating the 
time it takes for a particle to be mirrored or transmitted in the null electric 
field frame. In this frame the time rate of change of the particles parallel 
velocity V’ is given approximately by 


where V* is the particle gyration speed iA the null electric field frame, and 
1 

V* Is positive when the particle is directed towards the upstream region# 

H 

The effective V | B | experience by the particle at the shock discontinuity 
Is approximately the change in |b| divided by twice the particle gyroradlus 


V |B| ~ q | B | 

2mV' 


(26) 


Substituting equation (25) into (26) gives! 


dV' 

dt 11 


y.\ q (cos fk + ££!Lj!l) 

a 1 i 1 n 


where from the continuity of the normal component of B across the shock 

—1 9 9 

cos ij >2 ■ N cos In the null electric field frame V' + V' * 

H 1 

constant * V'2 which gives 


dV’ 


/v» 2-y« 2 


n (N-l) 

n«* sake* : it ..n cos fit (l + JL) dt 
8 1 N 


(27) 


where ^ «* q | f m"” 1 * . 

Integrating equation (27) gives 


n,(N-i) i 

a' ~ a' - — = cos \J»i (1 + i.) t 

18 1 N 


where a* » cos’" 1 V’/V’ and ot* is the initial value of the pitch angle in the 

II i 

null electric field frame. 

Using the fact that the pre- and post-shock values of V f are equal for 
reflected particles, and that the post reflection value of «’ » if - od , the 
shock interaction time for a reflected particle t^ is given approximately by 
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8 N sec 

2 (2a - it) . (2a) 

fy(N ~i ) i 


The relationship between a* ana the initial pitch angle in the upstream plasma 

i 

rest frame a is 
i 


cos «' * (cos a - R) (1 - 2R cos a + r2)-1/2 . 

1 i i 


Equation (28) above shows that tg is proportional to the particle’s mass 

to charge ratio rn/q. This result is independent of the assumption made in 

Equations (25) and (26). It comes from the fact that in the dimensionless form of 

the Lorentz force equation the dimensionless unit of time is proportional to 

ra/q. For a particle that is transmitted downstream the post transmission value 

of «’ * sin“^(* / N sin a 1 ) and the particle shock interaction time Un for such 

1 

particles is approximately 

8 N sec . 

trprv ~ * — — *■ *- [sin 1 (^N sin a') - a’] (29) 

TD n L (N 2 - i) i i 

For a particle that is transmitted upstream, the Sign of the right hand term in 

equation (25) is positive, the post transmission value of a* » sin"* (N^^sin/) 

and the shock Interaction time tqqj for such a particle is approximately 

8 N sec lji, -i /-> 

bmtf ~ 5— [sin 1 (N 1/z sin a') - a’] , (30) 

TtI %(N 2 - 1) 1 i 


Equations (28) ** (30) are not valid for shocks in which V^ x t an - 
+ * 

Jv s x n >| is greater than the speed of light, as there is no inertial Frame 
in which the shock induced E * 0. For this i|j 2 90* situation the interaction 
time can ho estimated from Equation (31), 

AT - J*E t ± (31) 
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where AT is the particle fractional change in kinetic energy * (N <* 1) 0.5 tnV z , 
is the grad B drt£t induced current|j| M qmV z 7|Bj(2q|B| z ” 1 ) t V^ is the particle 
gyration speed in its guiding center rest frame, and t j_ is the shock interaction 
time for a perpendicular shock. Combining the above terms gives, 

t, - !1" 1 (N+ 1)V, V‘ l (32). 

* 1 **“ lx . 
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V, NUMERICAL SIMULATION RESULTS 


A. Test of the Equality of Pre- and Post-Encounter p Assumption 


In order to test the assumption that the pre- and post-encounter values of 

p ara equal In the null electric field frame, It Is necessary to follow particles 

trajectories throughout the shock interaction in that frame and compare the 

ensemble average of the post-interaction value of p with the initial value. 

Since the particles' trajectory in the null electric field frame (and 

shock rest frame) can he expressed in terms of analytical functions, numerical 

integration techniques are not needed. The procedure is to choose the initial 

* # * 

position (x^, y_£, z-^) and velocity (x^, y^, »-<) of the particle and time step 


size At and then compute the position and velocity of the particle at time 


_ A 
-1 ” 


VUrt AH *%*%*’& 0*1 Art! « h j« | ma h -» rtwa t*lmn jiart/l no I t b it 1 {-• 

uic pv/SJ A V.-AVM* vVAVVA V/ l. imv v»*vu wwLv* uo vV/MUi 


tions to calculate the position and velocity at time 1*2 “ 2 At, and so on. 

The algorithms used to compute the particles' position and velocity are 
theoretically exact and in practice accurate to the single precision (li digit) 
accuracy of a Control Data 3800 computer on which the calculations were carried 
out. Positional and velocity errors occur if, during a step, the particle 
crosses the shock. The errors occur because gyroradius, gyrofrequency and 
drift velocity do not change during a step but do change across the shock. 

This error is minimized by using an iteration process to calculate the time 
the particles cross the shock and changing the step size accordingly so that 
the particle does not cross the shock but "ends up" at the shock surface and 
then starts the next step with the appropriate gyro radius, gyrofrequency, and 
drift velocity. 

That the pre- and post-encounter value of p for a particle which Is trans- 
mitted or reflected by a fast move shock wave should be equal is not self-evident 
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According to adiabatic theory of charged particle motion (Northrop, 1963], 
p averaged over a gyroperiod is conserved for static fields if 

L » Rg (33) 

where t is the scale length over which the magnitude of the magnetic field 
changes and R g is the particle gyroradius. Tn the shocks to be considered, ( B | 
changes discontinuously by a factor of 2 or more and clearly Equation (33) is 
violated. However, Pe3*es (1981] has shown that for perpendicular Bhocks 
('{'l ■ 90°), the equality of the pre- and post-encounter values of p is due to 
the continuity of the flux of particle plus field angular momentum through the 
shock. 

The equality of the pre- and post-encounter p assumption was tested by cal- 
culating Pfinal/hjLnitial for different values of the particle speed and pitch 
angle in the null electric field frame. 

The results for particles initially upstream for the case N * 2, t|^ * 89° 
and |V 1( » 10 Vi x are presented in Table 1. Column 1 presents the initial pitch 
angle where in this table, and this table only, , particles with a| < 90° 
are directed towards the shock. Column 2 gives the gyrophase averaged value 
of the ratio of the post-encounter value of p to the pre-interaction value 
( */ ¥)• Column 3 presents the gyrophase averaged fractional change in kinetic. 
Column 4 (5) gives the maximum (minimum) value of Pf/pj. for each pitch 
angle group. Column 6 gives tha range of the number of shock crossings. 

Column 7 tells whether particles within that pitch angle group were transmitted 
(T) or reflected (R). 

As E « 0 in the frame in which the simulation is carried out, the nonzero 
values of AT/Tj indicate that computer round off errors are occurring and/or 
the interation process to calculate the shock crossing time is not converging 
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fast enough. However, the errors are only l part in 10& even after 174 shock 

crossings and give confidence to the numerical technique used. Column 2 shows 

the precision with which the Initial and final gyrophase-averaged y are equal 

increases with increasing a', ranging from a difference of 8% at o’ * 10* to 

i i 

0.0008% at o^ - 81°. Columns 4 and 5 show the the precision with which an 

individual particle’s pre- and post-shock y are equal increases with increasing 

o’, with lower limits ranging from 68% at o' *■ l* to 0.01% at a' * 81*. Column 
i i i 

7 shows that test particles with o^ < 41* (> 51°) are transmitted (reflected). 
The latter result is consistent with the reflection/ transmission criteria 
discussed in section XV A2 which predicts that for a strength two shock upstream 
particles with ftj < (>) 45° are transmitted (reflected). 

The equality of initial and F inal y assumption for particles originating 
in the downstream, high }b| region Was » As in the initially up- 

stream case, the pre- and post-shock values of y were found to be equal to * \% 

for ct’ * 11° to .001% for a’ * 81°. Runs were also made for both upstream 
i i 

and downstream particles with Vj[ ** 20, 30 V^ x and with N « 4 with similar 
results. 


B. Comparison of Analytical Model Predictions of 
Energy Gains with Numerical Simulations 

1. Perpendicular Shocks 

Equation (14) predicts that in the limit ^ ( a perpendicular shock) 

AT/Tj^ * (N - 1) sin^ c^. This means that averaged over phase angle y f / * \. 
Numerical simulations of the interaction of charged particles with perpendicular 
shocks have been carried out previously. Parker [1958] found that averaging 
over entrance gyropbase the pre- and post-shock valus of " are equal for 
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infinitesimally thin fast-mode perpendicular shocks. Chen and Armstrong [ 1972 J 
and Posses l 1979), with more extensive studies, confirmed Parker's [1958] results. 
A typical particle trajectory in the x - y plane of a perpendicular shock is 
shown in figure Z of Pesses [1981]. 


2. Oblique Shocks 

Numerical studies of the Interaction of charged particles with oblique 
shocks have been done by Hudson [1965], Quenby and Webb [1973], Chen [1975], 
Pesses [1979], Terasawa [1979] and Decker [1981]. The first two papers are 
concerned with the conditions for reflection and transmission. Chen [1975] 
made a detailed study of the post-interaction energy and pitch angle distri- 
bution as a function of N, ^ and V^ X V^""^ . Chen also considered the effects 
of charged particle scattering by magnetic field irregularities. Terasawa 
[1979] did a study similar to Chen [1975] and also considered the effects of a 
finite shock thickness. 

The predictions of the analytical model of shock acceleration derived in 
Section IV have been compared with results from independent numerical simulations 
by Pesses [1979] which uses the same particle following technique described in 
Section \F A. 

The particle's initial position, velocity* pitch angle, and phase angle are 
Specified In a plasma rest frame and then transformed to the shock rest frame. 
After the particle-shock interaction is completed, the particle velocity vector 
is transformed back to the appropriate plasma rest frame. The results presented 
are averaged over initial and final gyrophase. 

In Figure 3, Equations (12) and (11) are compared to the numerical simu- 
lations results for the case ^ “ 84°, N “ 2 and Vi * 40 V-^. The vertical 
axis gives AT/T^, the horizontal axis o^. The portion of the curve labeled R 


- 25 - 


(T) [!.J is for reflected (transmitted) f noninteracting] particles. The analytical 
and numerical values are in excellent agreement. Passes [1979] has shown the 
agreement between the predictions of the analytical equations presented in 
Section IV of this paper and the numerical simulations for final pitch angles 
and reflection, transmission and noninteracting coefficients are alao excellent. 

G. Acceleration Time 

Equations (28) - (31) show that the particle's acceleration time in the 
shock is proportional to its gyroperlod* Hence the acceleration time is inversely 
proportional to the magnetic field strength and proportional to the particle 
mass to charge ratio m/q. 

To check the accuracy of the analytical expressions for acceleration time 

(Equations 28-31) the particle-shock interaction time is calculated numerically. 

The partial shock interact time is defined as the time between the particles 

first and last crossing of the shock front. The tJ"4C ns a function of the 

particle gyrofrequency in the upstream region JT^-is calculated using the same 

JL 

numerical procedure as described above for the Uf/bi* end AT/T-^ calculations. 

An example of the agreement between one numerical and analytical acceleration 
times is given in Figure 3. Tin. numerical (solid circles) and analytical cal- 
culated (solid curve) value of t^ and t r j) as a function of uj. are compared for 
the case; “ 82°, N * 2, and R * 1. The agreement between the analytically 
and numerically calculated acceleration time is very good for transmitted 
particles, and fair for reflected particles. 

The perpendicular shock interaction times pi*edictd by Equation (32) are 
nlso consistent with those calculated numerically by Passes [1979). For example, 

for N - 2 and V i V" 1 * 20 Pesses [1979] finds that t, « 42 ft." 1 , while Equation 
I- lx 

(32) predicts tj_ ~ 40 
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VI. CONCLUSION 


This paper has dealt primarily with the physics and calculations of the 
energy gain of energetic particles that are reflected or transmitted by mag- 
netosonlc fast-mode shock waves. For a typical interplanetary shock of magnetic 
strength N - 2 the maximum increase in a particle energy from one shock encounter 
is a factor of 5.83. For an ion with an initial energy of 30 keV/nuc to end up 
with 10 MeV/nuc at least four shock encounters are needed. Clearly, in order to 
understand interplanetary shock acceleration phenomena, which routinely result in 
10 MeV protons, it is necessary to extend the single encounter model in this 
paper to a multiple encounter model. Such a model is presented Paper II. 
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PRECEDING PAGE BLANK NOT FILMED 

FIGURE CAPTIONS 


Figure 1: Shock geometry In the shock rest frame. 


Figure 2: 

a and b: Graphic view of acceleration process at shock front in the upstream 

(a) and downstream (b) plasma rest frames. 

Figure 3: The dependence of A T/T^ on oi£, for <* 84° , N 3 2, = 40 V^ x . 

Model predictions shown as solid lines, numerical results shown as 
solid circles. R,T and N stands for reflected upstream, transmitted 
downstream and noninteracting particles, respectively. 

Figure 4: The dependence of acceleration Lime on for ■“ 82°, N = 2, and 

R ■ 1. Model predictions shown as solid lines, numerical simulation 
results shown as solid circles. t^d stands for reflected up- 

stream, transmitted downstream particles, respectively. 
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Table 1 


Conservation of First Adiabatic Invarient 




N - 2 

% - 89° 

* 89.5® 



a 

<Uf/Mi> 

at/t* 

Pf/ M^max 

Vif/p^min 

// Crossings 

R/t ' 

1 

1.08391 

1.096 x 10“ 6 

1.68348 

0.390064 

1 or 3 

T 

11 

1.00122 

1.366 x 10“ ? 

1.02527 

0.968477 

11 or 12 

T 

21 

0.999978 

2.565 x 10“ 7 

1.00765 

0.986034 

23 or 25 

T 

31 

0.999954 

3.738 x 10“ 7 

1.00233 

0.993105 

41 or 43 

T 

41 

0.999958 

5.431 x 10“7 

1.00191 

0.995753 

73 or 75 

T 

51 

0.999283 

9.238 x 1G“ 7 

1.00104 

0.997076 

174 

R 

61 

1.00005 

1.304 x 10" 6 

1.00196 

0.998086 

132 or 134 

R 

71 

1.00002 

1.248 x 10“6 

1.00059 

0.999567 

104 or 106 

R 

81 

0.999992 

1.199 x 10“6 

1.00009 

0.999891 

78 

R 
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